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Abstract. We present an universal treatment for a substantial nuclear suppression representing a 
common feature of all known reactions on nuclear targets (forward production of high-p^ hadrons, 
production of direct photons, the Drell-Yan process, heavy flavor production, etc.)- Such a sup- 
pression at large Feynman xf, corresponding to region of minimal light-cone momentum fraction 
variable %i in nuclei, is tempting to interpret as a manifestation of coherence or the Color Glass 
Condensate. We demonstrate, however, that it is actually a simple consequence of energy conserva- 
tion and takes place even at low energies, where no effects of coherence are possible. We analyze 
this common suppression mechanism for several processes performing model predictions in the 
Q} ' light-cone dipole approach. Our calculations agree with the data. 
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g '■ INTRODUCTION 



In the proton(deuteron)-nucleus and nucleus-nucleus collisions, investigated at the Rela- 
tivistic Heavy Ion Collider (RHIC), recent measurements of high-p^ particle spectra by 
the BRAHMS [1,2], STAR [3] and PHENIX [4] Collaborations show a strong nuclear 
suppression. Observed nuclear effects occur not only at large forward rapidities [1, 2, 3] 
but unexpectedly also at midrapidities [4] . 

Besides, quite strong and universal nuclear suppression at large Feynman xf is con- 
firmed by the collection of data from [5] for the production of different species of par- 
ticles in p — A collisions. The rise of the nuclear suppression with xf is also supported 
by the NA49 data [6] at lower energy corresponding to c.m.s. energy y/s = 17.3 GeV. 
The onset of strong nuclear effects at large pj has been also demonstrated for direct 
photon production in Au — Au collisions at RHIC by the PHENIX Collaboration [7]. 
The E772 experiment at Fermilab [8] first observed that the Drell-Yan (DY) process is 
considerably suppressed at large xf- 

Assuming large forward rapidities, the basic explanation for such an effect has been 
based on an idea that in this kinematic region corresponding to the beam fragmentation 
region at large Feynman xp one can reach the smallest values of the light-front momen- 
tum fraction variable X2 in nuclei. It allows to access the strongest coherence effects such 
as those associated with shadowing or the Color Glass Condensate (CGC). 

It was shown in refs. [9, 10, 1 1] that a considerable nuclear suppression for any reac- 
tion at large xf (small %2) is caused by another effects, which can be easily misinterpreted 



as coherence. Such a suppression can be treated, alternatively, as a Sudakov suppression, 
a consequence of a reduced survival probability for large rapidity gap (LRG) processes 
in nuclei, an enhanced resolution of higher Fock states by nuclei, or an effective energy 
loss that rises linearly with energy. It was demonstrated in refs. [9, 10] that the nuclear 
suppression at large xf is a leading twist effect, violating QCD factorization. 

In this paper we will analyze nuclear suppression at large rapidities (large xf) for the 
following processes occurring in p(d) — A and A —A collisions: 

• production of leading hadrons with small pj 

• high-px hadron production at forward rapidities in p(d)-A collisions 

• production of hadrons at small energies vs. NA49 data 

• high-pj hadron production at midrapidities 

• direct photon production in Au-Au collisions 

• Drell-Yan production at large xf 

SURVIVAL PROBABILITY OF LARGE RAPIDITY GAPS 

Treating any hard reaction, which is LRG process in the limits — » 1, gluon radiation is 
forbidden by energy conservation. If a large-.*^ particle is produced, the rapidity interval 
to be kept empty is Ay = — ln(l —xf). Assuming as usual an uncorrected Poisson 
distribution for gluons, the Sudakov suppression factor, i.e. the probability to have a 
rapidity gap Ay, becomes 

S(Ay) = e -< n G(Ay)> , (1) 

where nc(Ay) is the mean number of gluons that would be radiated within Ay if energy 
conservation were not an issue. 

The mean number (nc{Ay)) of gluons radiated in the rapidity interval Ay is related 
to the height of the plateau in the gluon spectrum, (Ho(Ay)) = Aydnc/dy. Then, the 
Sudakov factor acquires the simple form, 

S(x F ) = {\-x F ) dncldy . (2) 

The height of the gluon plateau was estimated in ref. [12] as, 

**l = *!lJA-\. (3) 

d y n \ a qcdJ 

For further calculations we take a s = 0.4 (see discussion in ref. [9]), which gives with 
high accuracy dno/dy = 1, i.e. the Sudakov factor, 

S(x f ) = 1-xf . (4) 

One can formulate nuclear suppression as xp — ► 1 as a survival probability of the 
LRG in multiple interactions with the nucleus. Every additional inelastic interaction 
contributes an extra suppression factor S(xp). The probability of an n-fold inelastic col- 
lision is related to the Glauber model coefficients via the Abramovsky-Gribov-Kancheli 



(AGK) cutting rules [13]. Then the survival probability at impact parameter b reads, 

A 



Wlbib) = & x V [-a!fT A (b)]t^\o^T A (b)} n S( XF ) n -\ (5) 



n=\ 



where T A (b) is the nuclear thickness function. 



PRODUCTION OF LEADING HADRONS WITH SMALL p T 

The left panel of Fig. 1 shows the collection of data from [5] for production of different 
species of particles in p —A collisions exhibiting quite a strong and universal suppression 
at large xp . Moreover, these data cover the laboratory energy range from 70 to 400 GeV 
and demonstrate so the xf scaling of nuclear effects. 

It is natural to relate the observed suppression to the dynamics discussed in the 
previous section. The nuclear effects can be calculated using Eq. (5) summing over the 
number of collisions and integrating over the impact parameter, 

R A/N (x F ) = \ fd 2 be- (Je ff TA ^ \e (l-XF)o eff T A (b)_ l X (6) 

1 (1 — XF)O e ffA J I J 

In the Glauber model G e ff = of n N . However, Gribov's inelastic shadowing corrections 
substantially reduce a e ff [14, 15]. 

To compare with data, the nuclear effects are parametrized as R a /n ^ A a , where the 
exponent a varies with A. We used A = 40, for which the Gribov corrections evaluated 
in [15] lead to a e ff ~ 20 mb. Then a simple expression Eq. (6) explains the observed xp 
scaling and describes rather well the data. 



HIGH-pt HADRON PRODUCTION AT FORWARD RAPIDITIES 

The cross section of hadron production in dA(pp) collisions is given by a convolution of 
the distribution function for the projectile valence quark with the quark scattering cross 
section and the fragmentation function, 



d 2 o v f, , f 2 d}a\qA{p^ 
-^^ = LJ dzf q/d{p) (x u q T ) -^—^ 

Zmin 



d2 ° v- f j r / 2 ^d 2 o[qA(p)\ 

D h/q{Z), (7) 

q ^ " ~u "-i qr=PT/z 



where x\ = ^e* 1 . The quark distribution functions in the nucleon have the form using 
the lowest order parametrization of Gluck, Reya and Vogt [16]. We used proper frag- 
mentation functions using parametrization from [17]. 

Interaction with a nuclear target does not obey factorization, since the effective projec- 
tile quark distribution correlates with the target. The main source of suppression at large 
Pt concerns to multiple quark rescatterings in nuclear matter. Summed over multiple 
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FIGURE 1. (Left) Exponent describing the A dependence (°= A a ) of the nucleus-to-proton ratio for 
production of different hadrons as a function of xp. (Right) Ratio of negative hadron and neutral pion 
production rates md—Au and pp collisions as function of pj at pseudorapidity 77 = 3.2 and 77 = 4.0 vs. 
data from the BRAHMS [1] and STAR Collaborations [3], respectively. 



interactions, the quark distribution in the nucleus reads, 



f q / N {x\AT) = Cf q /N(x\,q 2 T ) 



fd 2 b 



e -x\O ef fT A (b) _ e -(J ef fT A (b) 



(\- Xl )Jd 2 b 



1 _ e -o e ffT A {b) 



(8) 



where the effective cross section o e ff = o e ff{pT,s) = {o^ q (rr)) / {Oq q {r T )) has been 
evaluated in [9]. The normalization factor C in Eq. (8) is fixed by the Gottfried sum rule. 

The cross section of quark scattering on the target do[qA(p)]/d 2 qTdr\ in Eq. (6) is 
calculated in the light-cone dipole approach [18, 19]. In our calculations, we separate 
the contributions characterized by different initial transverse momenta and sum over 
different mechanisms of high-pr production. Details can be found in [9]. 

The BRAHMS Collaborations [1] in 2004 found a substantial nuclear suppression 
for high-/?7 negative hadrons produced at pseudorapidity r\ = 3.2. Two years later, the 
STAR Collaboration [3] has been observed even stronger suppression for neutral pions 
at r\ — 4.0 as one see from the right panel of Fig. 1 . Because the data cover rather small 
%2 ~ 10~ 3 , the interpretation of such a suppression has been tempted to be as a result of 
saturation [20, 21] or the CGC [22], expected in some models [23]. 

Even if one supposes to interpret the observed suppression at 77 = 3.2 in terms of 
CGC, such an interpretation should fail at larger r\ = 4.0, where the observed suppres- 
sion is more than a factor of 2 larger. The stronger onset of the quantum coherence 
effects at r\ = 4.0 can not explain such a huge rise of nuclear suppression. 

Much stronger nuclear effects at r\ = 4 can be simply explained by the energy 
conservation as a much smaller survival probability of LRG at larger T]-values [9, 1 1]. 

Energy conservation applied for multiple parton rescatterings leads to xp scaling 
of nuclear effects [9, 10, 11]. We expect approximately the same nuclear effects at 
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FIGURE 2. (Left) Theoretical predictions for an approximate exp(rj)/^/s- scaling of the ratio 
Rd+Au(Pr) for n° production rates in d—Au and pp collisions. (Right) Ratio, R p+ pi,(pt), for 7T* pro- 
duction rates in p — Pb and pp collisions as function of pj at two fixed values of Feynman xf = 0.025 
and 0.375 vs. the NA49 data [6]. 



different energies and pseudorapidities corresponding to the same values of xp. Such 
a situation is demonstrated in the left panel of Fig. 2, where we present pp dependence 
of nuclear attenuation factor R-d+AuiPr) f° r K° production at different c.m.s. energies 
and r\ keeping the same value of xp . 

NUCLEAR SUPPRESSION AT SMALL ENERGY VS. NA49 DATA 

The right panel of Fig. 2 clearly demonstrates a stronger onset of nuclear effects at larger 
xf. The model predictions for nuclear suppression have been performed employing the 
dipole formalism and using the mechanisms for the valence quarks described in [9]. One 
can see a good agreement of our calculations with NA49 data [6] . 



mGH-p T HADRON PRODUCTION AT MIDRAPIDITIES 

As a consequence of xp- scaling is an expectation of similar nuclear effects also at 
midrapidities. However, the corresponding values of pj should be high enough to 
keep the same value of xp. Such an expectation is confirmed by the recent data from 
the PHENIX Collaboration [4] showing an evidence for nuclear suppression at large 
pj > 8 GeV (see the left panel of Fig. 3). 

At r\ = the small-p^ region is dominated by production and fragmentation of 
gluons. On the other hand, the region of very large pp is dominated by production 
and fragmentation of valence quarks. Consequently, any value of the hadron transverse 
momentum differs only in the relative contributions of valence quarks and gluons. 

It means that we include also gluons in our calculations. Details can be found in 
ref. [24]. Correspondingly, the cross section for hadron production, Eq. (6), is extended 
also for gluons with corresponding distribution function, parton scattering cross section 



and the fragmentation function. Including multiple parton interactions, the gluon distri- 
bution in the nucleus is given by the same formula as for quarks (see Eq. (8), except 
o e ff, which should be higher by the color factor 9/4. 

If the effects of multiple parton rescatterings are not taken into account the pr 
dependence of R(1+Au{pt) is described by the thin dashed line. One can see from the 
left panel of Fig. 3 that our calculations at moderate pr are not in a bad agreement with 
data and a small suppression at large pj is given by the isospin effects. After inclusion 
of multiple parton rescatterings the model predictions presented by the thin solid line 
underestimate the data at moderate pj. However, at larger pj quite a strong onset of 
nuclear effects is not in disagreement with corresponding experimental points. 
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FIGURE 3. (Left) Nuclear attenuation factor Rd+Au(pr) as a function of pj for production of n° 
mesons at ^fs = 200 GeV and T] = vs. data from PHENIX Collaboration [4]. (Right) Nuclear modi- 
fication factor for direct photon production in Am —Au collisions as a function of p j. 

Calculations in the RHIC energy range at midrapidities are most complicated since 
this is the transition region between the regimes of long (small pj) and short (large 
Pt) coherence lengths. Instead of too complicated rigorous light-cone Green function 
formalism [25, 26, 27, 28] we preset corrections for finite coherence length using the 
linear interpolation performed by means of the so-called nuclear longitudinal form 
factor [24]. Such a situation is described by the thick solid and dashed lines reflecting 
the cases with and without inclusion of the multiple parton rescatterings, respectively. 
It brings the model predictions to a better agreement with data at moderate pj. Nuclear 
suppression at large pj > 10 GeV observed by the PHENIX experiment [4] can not be 
explained as a result of CGC because data cover rather large X2 ~ 0.05 — 0.1. 

DIRECT PHOTON PRODUCTION IN AU-AU COLLISIONS 



Expressions for the production cross sections have been derived employing the dipole 
formalism [29, 30, 31, 19, 32]. Model predictions for Rau-Au as a function of pr are 
compared with the PHENIX data [7] in the right panel of Fig. 3. If multiple parton 
rescatterings are not taken into account the model calculations depicted by the dash- 
dotted line overestimate the data at large pj J> 13 GeV. The onset of isospin effects gives 



a value Rau-Au —* 0.8 in accord with our calculations. Inclusion of the multiple parton 
rescatterings leads to a stronger nuclear effects at large pj as is demonstrated by the 
dashed line. It brings a better agreement of the model with data. Finally, the solid line 
additionally includes also a small correction for the EMC effect [33]. 

DRELL-YAN PRODUCTION AT LARGE x F 

The DY reaction is also known to be considerably suppressed at large xf [34] as one 
can see from Fig. 4. Model calculations have been performed using expressions for the 
production cross sections in the color dipole approach [30, 31]. We included also the 
effect of multiple parton rescatterings [9, 10, 11] discussed above. Model predictions 
are in a reasonable agreement with data from the E772 experiment [8]. 





FIGURE 4. Ratio of Drell-Yan cross sections on Tungsten and Deuterium as a function of x\ 



SUMMARY AND CONCLUSIONS 



In this paper we analyze a significant nuclear suppression at forward rapidities (large 
xf) for several processes. The new results are the following : 

• QCD factorization fails at the kinematic limits, xp — ► 1, x\ — ► 1. Nuclear targets 
cause a suppression of partons with x — ► 1, due to energy sharing problems. 

• Suppression of high-pr hadrons at large rapidity observed by the BRAHMS and 
STAR Collaborations is well explained. 

• We predict x\ (xf) scaling, i.e. the same nuclear effects at different energies and 
rapidities corresponding to the same value of x\ {xf). 

• Model predictions are in a good agreement with NA49 data [6] and clearly demon- 
strate the rise of nuclear suppression with xp. 

• Predicted strong nuclear suppression for the large-p^ direct photon production in 
Au —Au collisions is in a good agreement with the PHENIX data [7]. 

• According to xp scaling we predict nuclear suppression at large pj also for hadron 
production at r\ = 0. Model calculations describe well the PHENIX data [4]. 



• Study of nuclear effects at midrapidities is very important because at large pj the 
data cover rather large %2 ~ 0.05 — 0.1, where no effect of coherence is possible. It 
allows to exclude the saturation models or the models based on CGC. 

• Suppression of Drell-Yan pairs at large xp observed by E772 Collaboration [8] is 
well explained. 
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